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In a recent letter, we have studied the scattering of relativistic electrons (or positrons) by 
the Coulomb field of a nucleus in two dimensions [1]. The Dirac equation was solved in polar 
coordinates, and the scattering amplitude was obtained as a partial wave series. For light 
nuclei the series can be summed up approximately to give a closed result. The result, though 
being approximate, exhibits some nonperturbative features and cannot be obtained from the 
lowest-order contribution of perturbative quantum electrodynamics (QED). This feature is not 
manifest in the corresponding result in three dimensions. The purpose of the present paper 
is to extend our previous work to the case with both Aharonov-Bohm (AB) and Coulomb 
potentials. In the nonrelativistic case, this can be described by the stationary Schrodinger 
equation 



where (r, 6) are polar coordinates on the xy plane, and has been studied in the literature [2, 
3]. Note that this is a two-dimensional model and is different from the three-demensional 
Aharonov-Bohm-Coulomb (ABC) system studied in the literature [4-9]. Of course the AB 
potential is the same in either two or three dimensions. The difference lies in the Coulomb 
field. Previously the ABC system (in two or three dimensions) is regarded as describing a 
charged particle (with charge —q) moving in both an AB (with fiux 0) and a Coulomb field 
generated by external sources. Thus if the Coulomb field is generated by a nucleus (with 
charge Zq), one should have a magnetic fiux string fixed on the nucleus to generate the AB 
potential. The situation seems difficult to realized in practice. In a recent paper [3], we 
have shown that the above equation can describe the relative motion of two particles, one 
carrying electric charge and magnetic flux (— g, —(p/Z) and the other {Zq, 0), where Z 
is a real number. Then n is the reduced mass of the system. Such particles appear in (2-1-1)- 
dimensional Chern-Simons field theories as charged vortex soliton solutions [10-17]. They are 
also called anyons because their angular momentum may take on values other than integers or 
half integers. Applications of such objects can be found in the study of fractional quantum Hall 
effect [18-21], superconductivity [21-22], and repulsive Bose gases [23]. Of course the solitons 
have finite sizes, and the real interaction between them may be rather complicated. Thus 
the above description by the two-dimensional ABC system is merely a rough approximation. 
Anyway, the model may be of interest in itself since exact analysis is possible. A relativistic 
case described by the Klein-Gordon equation has been studied in the literature [24]. Here we 
will deal with a spin ^ particle described by the Dirac equation. 

Let us begin with the stationary Dirac equation translated from Eq. (1): 



where k. = Zq"^, p is the momentum operator, a = 7^7, and 7^ = (7^,7) {n = 0, 1,2) are 
Dirac matrices satisfying {7'^,7''} = 2g^^'^ where g^'^ — diag(l, —1, —1). In two dimensions the 
Dirac matrices can be realized by 2 x 2 matrices: 



where the a's are Pauli matrices. Thus ip is a two-component spinorial wave function. As in 
the case with a pure Coulomb field, it is not difficult to show that the conserved total angular 
momentum operator is 




(1) 





7° = cr^, 7^ = icr\ 7^ = icr^. 



(3) 
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J = X'Py- yp:, + y7S^- (4) 

Therefore the particle has spin |. 

The Dirac equation (2) can be solved in the polar coordinates by separation of variables. 
Bound-state solutions (for Z > or k, > 0) are relatively easy to obtained, and some results 
will be given at the end of the paper. Scattering solutions exist when E > fic? oi E < —fic^. 
The latter correspond to antiparticles after second quantization. At the level of single-particle 
theory, their scattering can be treated formally in a way similar to the former. Thus we only 
consider scattering solutions with E > jic^ in the following. We use the representation (3). 
Let 

mo + 1', (5) 



27r/ic 

where mo is an integer and — 5 < < 5, and 



'^/(r)exp[i(i-mo-l/2)^]/V^^ 
y g{r) exp[i(j - mo 1/2)6*]/ ^/r j 



4,4.... (6) 



It is easy to show that Jijjj = {j — mo)?iijjj. Thus j is a good quantum number. Substituting 
this expression into Eq. (2) we obtain two coupled ordinary differential equations for the radial 
wave functions: 

/ + kig + = 0, (7a) 

ar r r 

— H g - k2f - -/ = 0, 7b 

ar r r 



where 



^^^EjU^ ^^^^JZJ^ ^=J1^^ (8) 
he ^ he ^ he he 



Then we introduce the new variable 



he 

and two new functions u{p), v{p) through the definition 



(9) 



fir) = ^^/he'nu{p)+v{p)], g{r) = -'-^,e^P[u{p) - v{p)], (10) 



to recast Eq. (7) into the form 



du ip ip' + j + u 

u v — 0, (11a) 

dp p p 

dv id i3' — j — 1/ / , s 

— + 2iv + —v + — u = 0, (lib 

dp p p 



3 



where 

(12a) 




(12b) 

where Vc is the classical velocity of the incident particle. This is more convenient. Indeed, one 
can eliminate v immediately to obtain an equation for u alone: 

where we have used — = 7^. In this paper we assume for convenience that I7I < | (for 
electron- nucleus interaction, this means that Z < 68, which is in general satisfied in practice). 
This is not quite enough. If is not close to |, we further assume that I7I < || ± 
(This cannot hold if \v\ is very close to |, which causes some difficulty and will be discussed 
separately in the following. A larger 7 causes more difficulty.) Then for any j, the solution of 
Eq. (13) is well behaved at the origin. Let 



u{p) = p'w{p), s=^{j + uy--i^. (14) 
Then we have for w the equation 

(P'w , , dw , „ . / s 

2s + l + 2ip)— + 2(/3 + is)w = 0. (15 
op^ dp 

This is familiar. The solution that is well behaved at the origin is w{p) = ^{s—i(3, 2s-|-l, —2ip), 
where $(a, 6, 2;) is the confluent hypergeometric function [25]. So we have 

Uj(p) = ajp^^(s - i(3, 2s + 1, -2ip), (16a) 

where aj is a constant, and the subscript j of Uj that is omitted above has been recovered. 
Substituting this solution into Eq. (11a) we have 

vAp) = . ' ~'f-a, P'Hs + l-iP,2s + 1, -2ip), (16b) 
where we have used the formula 

+ $(a, b, z) = a$(a + 1, b, z), 

which can be obtained from other relations given in mathematical handbooks [25] . It should be 
remarked that ipj may be slightly singular at the origin when j' = ±|. However, the integral of 
ipjipj over any finite volumn converges and the solution is acceptable. This situation is similar 
to the case of a pure Coulomb field [1]. We take 

A^,/ ■ -^/x r(s — i/9) //Stt . .sn .7r\ , 

a,- = A2*(? + i/ + z/?')H ^exp — + zmvr - z hi- , (17a) 

^ ^ ^r(2s + 1) ^ V 2 2 4y ' ^ ^ 
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where m — j — ^ and 



A = u 



1 



then when r ^ oo we have for the radial wave functions the asymptotic forms 



(17b) 



cos i kr + pin 2kr 



mn 71 \ 1 



+ - {Sj — 1) exp ( ikr + ip In 2kr — i 



TT 



/ — r / ITITT TT \ 1 / TT \ 

5'j(^) ~^ ^ V ^2 ^™ sin Ikr + /Sin 2A;r ^ 4 / ~'~ 2^^'^-' ~ ^"^^ V ~'~ ^'^ ~ ^ 4 ) 

to the lowest order, where 



18a) 

)] 

(18b) 



Y(s - i(3) 

Sj = exp{2ir]j) = ( j + i/ + ^ ^ .^-^ exp(ij7r - isn) (19) 

and the rjj^s are phase shifts. The asymptotic form for ip — J^ji^j^ where the summation is 
taken over all j, turns out to be 



where 



with 



and 



+ 00 

E-r 
^ 

m=—oo 



exp(— imofc' 
\/ki + k2 



i\fk[ 



\ 



nkr 



cos ikr + (3ln 2kr — 



rriTT TT 



'^exp{tkr + if3ln2kr)f{e) """P^ '"""^^ 
r ^ Vki + k2 



- - e 



(20) 



(21a) 



(21b) 



:22] 



with 




X:exp(.r;,) sinr^.e^^^ = ^2(8, - l)e""^ 



(23) 



In spite of the presence of the AB potential, the probability current density associated with a 
solution ip is still given by 



j — cip^cxip. 



(24) 



As in the case of a pure Coulomb field [1] , one can show that the above t/j represents a correct 
scattering solution of which ipi^ is an incident wave and ipse the scattered one, f{9) is the 
scattering amplitude, and the differential cross section is given by 



a{9) = \f{9)\'. 



(25) 
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The choice of aj in Eq. (17) is thereby proved to be appropriate. Since J2j e*"*^ = 27r5{9), we 
have for 9 ^ 



f{0) 



Am6 



(26) 



This result with 5'^ given by Eq. (19) is exact. However, it is even more difficult to sum up 

the above partial wave series than in the case of a pure Coulomb field. 

Let us consider the case with 7^ <^ 1 and try to work out a closed result. For electron- 
nucleus interaction wc have 7 ~ Z/137, so the above condition means that Z is small, say, 
Z < 5. In this case we may approximately replace s by \j + i/| in Eq. (19). Note that (3 also 
depends on 7 and we do not make approximation with it. Thus the result will possess some 
nonperturbative features in regard to the parameter 7. With the above approximation, Sj is 
replaced by 

.r(m + i/ + l/2-i/3) _.^^r{m + u + 1/2 -i/3) 



r{m + u + l/2 + ip) 
{\m\-u + l/2-ip) 



t{(3 - (3')e- 



+ i{(3 - 



r(m + u + 3/2 + i(3)' 
.r{\m\ -u-l/2-ip) 



{j > 0), (27a) 



(j < 0). (27b) 



r(|m| -i/ + l/2 + i/3) ' r(|m| -i/ + l/2 + i/3)' 

It can be shown that the first term in either equation equals exp(2i5^) where 5^ is the non- 
relativistic phase shift of the mth partial wave, except for m = when — | < < 0. (The 
physical reason for the latter disagreement is not quite clear to us.) The second term is a 
relativistic correction which vanishes when Vc/c — > 0. Substituting Eq. (27) into Eq. (26), we 
have an approximate result for f{9): 



(28) 



where 

/o(^,^) 



V27rk 



.^^ r{i/2 + u-ip) 
r(i/2 + 1/ + i/3) 

T(3/2 -ly-ip). 



F{1, 1/2 + i/ - zP, 1/2 + u + 1(3, e'') 
e -F(l, 3/2 - i/ - 3/2- v + e"^^) 



V{3/2-v + i0) 

is the nonrelativistic partial wave result (somewhat different when — | < i/ < 0) and 



(29) 



h{9,v) 



V27rk 



-r(i/2 



t(3) 



-e" 



r(3/2 + 1/ + i/3) 

.T{l/2-u-i(3) 



F(l, 1/2 + u- ip, 3/2 + u + ip, e'") 
e '^F(l, l/2-u-ip,3/2-iy + ip, e" 



r{3/2-U + ipy - - ' -r^,- yj (30) 

is the relativistic correction. In these equations F{ai, 02, 61, z) is the hyper geometric function 
[25]. Using their functional relations it can be shown that 

. _,,g r(l/2 -u + tP)r{l/2 + v-iP) exp(^/jlnsin^ 9/2) 

m)v{i/2+zp) 

T{3/2-v -ip) _ 
T{3/2-v + ip) ~ 



^/2^^ 



'2ks\n9/2 
.T{~l/2 + v-ip) 
V{-l/2 + v + ip) 



xe''^F{l, 3/2- V- ip, 3/2- v + ip, e''^), 



(31) 
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and 



_ / _ ^\ r(l/2 + u- ip)T{l/2 -u + ip) e-'^/^-'"^ exp(z/3 In sin^ 9/2) 



V2nk 



.r(i/2 



r(3/2 



r(-l/2 + i/-z/3) 



r(i/2 



xe 



1/2 - z/ - il3, 3/2 - z/ + e"*') 



(32) 



where < ^ < 2tt. Note that these approximate expressions are well defined for the whole 
range of z/, and will be employed in the following discussions for the case when |z/| is close to 
i. Unlike the case with a pure Coulomb field, here we have additional terms in both /o(^, v) 
and fi{9, u), involving the hypergeometric functions. Thus closed forms are possible only for 
special z/'s when the additional terms vanish. This happens when z/ = and = |. Since the 
above discussions are not valid for |z/| close to |, we have now a closed result only when z/ = 0, 
or q(j)/27rhc takes on integer values. The result is the same as in the case of a pure Coulomb 
field: 



f{9) = r{9,0) = -t 



r(l/2 - i(3) exp(z/3 lnsin2^/2) 



2fcsin^/2 



(33) 



Since the result is more singular than 6{9) when 9^0, the S{9) term that is dropped above for 
9^0 can indeed be dropped everywhere and the above expression is enough. The differential 
cross section reads 



ai9) 



P tanh Ptt 
2ksm^ 9/2 



1 



sm 



9\ Ktanhlnn/hvc 



2/iw2 sin^ 9/2 



■ sm 




(34) 



where the first factor in the last expression is the exact nonrelativistic result and the subse- 
quent ones are due to the relativistic effect. Some remarks similar to those made in Ref. [1]: 
First, we have not made any approximation in regard to the incident velocity, so the result 
is valid for high energy collision. It is obvious that the relativistic correction becomes signif- 
icant when Vc is comparable with c. Second, though the above result holds for small 7 only, 
it involves a nonperturbative factor tanh/Jvr (note that /5 cx 7). Thus the result cannot be 
obtained from perturbative QED at the tree level. Moreover, if mo or (p is large, perturbative 
QED seems not applicable, but the above calculations hold as well. 



In the nonrelativistic case when z/ = and mo 7^ 0, we got an interference term in the cross section 
in additional to the result for a pure Coulomb field, because we excluded the s-wave solution which 
is slightly singular at the origin [3]. Now it seems that the s-wave is acceptable and that interference 
term is not necessary, because the potentials themselves are rather singular at the origin. Indeed, 
in the relativistic case, the solutions with j = ±1/2 are much more singular, and the singularity is 
essentially the same as that in a pure Coulomb field. 
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If y is close to | (— |) but I7I < |i ± z/|, then the above discussions are still valid. But then 
the replacement of S-1/2 {S1/2) by S^ii2 ('S'j^/2) is a poor approximation. In this case some 
corrections are necessary. Since no closed result is available, we do not discuss it in detail. 

Now we turn to the case when u is very close to | such that I7I > || — (the typical 
case is v = i), then the above discussions have to be modified. We will discuss it in some 
detail. (The other case with v very close to — | can be discussed in a similar way and will be 
omitted.) The crucial point is that s becomes imaginary when j = —\- (A larger 7 causes the 
same difficulty for some more values of j.) So the solutions of Eq. (13) are now given by 

u%2{p) = a%exp(27'lnp)<I>(-2/? + ^7', 1 + 2^7', -2ip), 



u 



(2) ^.^_.(2) 
-1/2 



(p) = aL^l/2exp(-z7'lnp)$(-z/5 - zf, 1 - 2i-f\ -2tp), (35) 



where 



7' = - (1/2 - i^r. (36) 



Here the two solutions have essentially the same singularity at the origin. They oscillate very 
rapidly and tend to no limit when p ^ 0. In the conventional opinion of quantum mechanics, 
such solutions are not acceptable [26, 27]. One possible way of resolving the problem is to 
cutoff the singular potentials (which represent some idealization) at some small radius. But 
this causes some mathematical difficulties. On the other hand, some authors have attempted 
to handle such solutions in some convenient way, and for bound states some proposals have 
been put forward to choose the appropriate linear combination of the two solutions [28, 29]. 
For scattering solutions the situation seems more involved. Here we will resolve the problem 
by some simple consideration. The solutions f i^^/2(p) "^-1/2 (p) '^^^ be easily obtained from 
Eq. (11a) and the above m's. We will not write them down. We take 



a^'^ -A 
"-1/2 — ^ 



p + ^ I' ] — -, ^ exp h ^7 m 2 — z— , 37 

^ \2 J\ r(l + 2i7') ^ \, 2 2 ' 4; ' ^ ' 



then the asymptotic forms for /l\''/2(^) ^^"^ 5'i^i/2(^) '^'^'^ given by Eq. (18) with 5'_i/2 replaced 

by 



-exp(7 7r). (38) 



r(i + i/5 + i7') 

-1/2 



However, 7^ 1, so S^^^ cannot be expressed as a phase factor. We then take af^]y2 by 

replacing 7' by —7' in a^]/2' tben the asymptotic forms for j^j^if) 5'i^i/2('") given by 
Eq. (18) with S'_i/2 replaced by S'^^ which is obtained by replacing 7' by —7' in S^^^ . As the 
two solutions above are equally preferable, we take the mean as the required solution. Then 
the asymptotic forms for /_i/2(r) and gi_i/2(r) are given by Eq. (18) with 5'_i/2 replaced by 

5= i[5W +5(2)]. (39) 

The solutions with j ^ are all the same as those given before. So we finally obtain the 
scattering amphtude (for ^ 7^ 0) 
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/(^) = --!=[ E S,e^-' + Se-% (40) 
V27r/c ^._,_i/2 

This is a very complicated result. When 7^ ^ 1, one can make approximations as before, but 
closed results are available only when v — \, or q(l)/27rhc takes on half integer values. We have 
then 

m = 1/2) - -^{S - St,/,)e-''. (41) 

Note that is well defined in Eq. (27b). Since we have neglected 7^ in calculating 

f^{9, 1/2), wc can as well neglect it in calculating S. It then turns out that S — Sti/2 to the 
first order in 7. So we have 



f(ft\ r-(fi 1 ,-ie/2 ^r(-i^) exp(i/jlnsin 9/2) 

m = f {6, 1/2) = -e r(l/2 + z/?) V2^sin^/2 

The differential cross section reads 



(42) 




Bcoth.371 / v"^ . o 0\ «;coth(7rfi;//ifc, , . 7 - i i -r i /.^\ 

= hn4r2 - # 2 j = 2K lin4/2 (' ' # 2 j (^ " # j ' 

where the first factor in the last expression is the exact nonrelativistic result [3], and the 
subsequent ones are due to the relativistic effect. The remarks made under Eq. (34) are also 
applicable here. 

Finally we give some results about the bound states. If I7I < || ± the solutions can be 
obtained without difficulty. The energy levels are 

= [l + 72/(n + s)2]l/2' 

where n = 0, 1, 2, ... is a radial quantum number, and s is defined in Eq. (14), which depends 
on j. The level with n = is not degenerate regardless of since one can show that only 
solutions with j > is possible in this case. When n > 0, the degeneracy depends on u. If 
z/ 7^ 0, there is no degeneracy. If z/ = 0, the energy level depends on |j| rather than j, and 
solutions with both positive and negative j exist. So the level is double degenerate. It is 
remarkable that there are no negative energy levels. The wave functions are given in terms of 
confiuent hypergeometric functions. Since the results are complicated we will not write them 
down. If|7| > ||±i/|, s will become imaginary for some j. For such values of j, special 
treatment [28, 29] of the solution is necessary, and the results are rather involved. We will not 
go into further details here. 

In conclusion, we have calculated the scattering amplitude and differential cross section 
for fast particles with ABC interaction in two dimensions. Exact results are given in partial 
wave series in general cases. Approximate results in closed forms are given in special cases. 
Though being approximate, the results exhibit some nonperturbative features and cannot be 
obtained from the lowest-order contribution of perturbative QED. We have also discussed the 
bound-state solutions and given some results. 

The author is grateful to Professor Guang-jiong Ni for encouragement. This work was 
supported by the National Natural Science Foundation of China. 
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